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We derive the spin state of the nucleus of Comet 103P/Hartley 2, its orientation in space, and its short- 
term temporal evolution from a mixture of observations taken from the DIXI (Deep Impact Extended 
Investigation) spacecraft and radar observations. The nucleus is found to spin in an excited long-axis 
mode (LAM) with its rotational angular momentum per unit mass, M, and rotational energy per unit 
mass, E, slowly decreasing while the degree of excitation in the spin increases through perihelion pas- 
sage. M is directed toward (RA, Dec; J2000) = 8 ± 4°, 54 ± 1° (obliquity = 48 ± 1°). This direction is likely 
changing, but the change is probably <6° on the sky over the ~81.6 days of the DIXI encounter. The mag- 
nitudes of M and E at closest approach (JD 2455505.0831866 2011-11-04 13:59:47.310) are 
30.0 ±0.2 m 2 /s and (1.56 ±0.02) x 10 _3 m 2 /s 2 respectively. The period of rotation about the instanta- 
neous spin vector, which points in the direction (RA, Dec; J2000) = 300 ± 3.2°, 67 ± 1.3° at the time of clos- 
est approach, was 14.1 ± 0.3 h. The instantaneous spin vector circulates around M, inclined at an average 
angle of 33.2 ± 1.3°, with an average period of 18.40 ± 0.13 h at the time of closest approach. The period of 
roll around the principal axis of minimum inertia (“long” axis) at that time is 26.72 ± 0.06 h. The long axis 
is inclined to M by ~81.2 ± 0.6° on average, slowly decreasing through encounter. We infer that there is a 
periodic nodding motion of the long axis with half the roll period, i.e., 13.36 ± 0.03 h, with amplitude of 
~1° again decreasing through encounter. The periodic variability in the circulation and roll rates during a 
cycle was at the 2% and 10-14% level respectively. 

During the encounter there was a secular lengthening of the circulation period of the long axis by 
1.3 ± 0.2 min/d, in agreement with ground-based estimates, while the period of roll around the long axis 
changed by ~-4.4 min/d at perihelion. M decreased at a rate of -0.038 (m 2 /s) per day in a roughly linear 
fashion. Assuming a bulk density between 230-300 kg/m 3 and a total volume for the nucleus of 
8.09 x 10 8 m 3 , the net torque acting on the nucleus was in the range 0.8- 1.1 x 10 5 kg m 2 /s 2 . In order 
to bring the spacecraft photometric and imaging data into alignment on the direction of M, the directions 
of the intermediate and short principal axes of inertia had to be adjusted by 33° (on the sky) from the 
values indicated by the shape model with an assumed homogeneous interior. The adjusted direction of 
the intermediate axis is RA, Dec = 302°, -16.5°. The morning and evening terminators in the images 
are identified, and the variation of the insolation at three regions on the nucleus associated with active 
areas calculated. The plume of water vapor observed in the inner coma is found to be directed close to 
the direction of local gravity over the sub-solar region for a range of reasonable bulk densities. The plume 
does not follow the projected normal to the surface at the sub-solar point. 

© 2012 Elsevier Inc. All rights reserved. 


1. Introduction 

* Corresponding author. Fax: +i 520 795 6286. Prior t0 the EP0XI ("DIXI”) mission encounter with the nucleus 

E-mail address: mbeiton@dakotacom.net (m.j.s. Belton). of Comet 103P/Hartley 2 (called 103P hereafter), near-aphelion 
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observations by Meech et al. (2009) indicated the presence of a 
periodicity near 16.6 h in the comet’s light-curve, and Lisse et al. 
(2009) had pointed out that the small size of the nucleus of the co- 
met (~0.57 km effective radius) combined with its relatively large 
perihelion water production rate (3 x 10 28 mol/s, A’Hearn et al., 
1995) implied that the spin state of the nucleus was likely highly 
excited and might substantially change in a single perihelion 
passage. 

Around the time of encounter (on 2011-11-04 with closest ap- 
proach at 13:59:47.310, i.e., JD 2455505.0831866), numerous 
Earth-based time-series observations of structures in the comet’s 
coma (Meech et al., 2011; Knight and Schleicher, 2010, 2011; 
Samarasinha et al., 2010, 2011; Jehin et al, 2010), radar observa- 
tions of the nucleus (Harmon et al., 2010, 2011), and high-spectral 
resolution microwave line observations of HCN (Drahus et al., 
2011) indicated that the magnitude of the periodicity was increas- 
ing with time and that the spin of 103P’s nucleus was probably in 
an excited state. That this was definitely the case was indepen- 
dently determined in the preliminary approach photometry and 
imaging results of the DIXI mission (A’Hearn et al., 2011; Belton 
et al., 2011a). While all of these studies were in broad agreement 
as to the magnitude and changes in the dominant periodicity and 
that a second periodicity was present, discussions at a special 
workshop on the subject held at the Kitt Peak National Observatory 
on April 5-6, 2011, showed less unanimity on the space direction 
of M (the rotational angular momentum per unit mass vector), 
the magnitude of the roll or oscillation periodicity, or whether 
the excited motion was in the short (SAM) or long (LAM) axis mode 
(see the Appendix in Samarasinha and A’Hearn (1991) for a de- 
tailed explanation of torque-free excited spin and the SAM and 
LAM modes). 

In this work we present a detailed analysis of DIXI observations 
that bear on the nucleus spin and show how the spin state changes 
through the DIXI encounter and that, while the rotation state is 
definitely excited, it is much less so than anticipated by Lisse 


et al. (2009). We also determine the average space direction of M 
and the motion of the instantaneous spin vector throughout the 
encounter. 

In Section 2 we collect all of the information about the nucleus, 
derived mainly from DIXI observations, that is needed for the anal- 
ysis. In Section 3 we present our basic assumptions and an over- 
view of the technique used to determine the spin state and its 
evolution. In Section 4 we analyze the approach and recession pho- 
tometry to determine the dominant periodicities in the light-curve 
and show that there are two groups of periodicities one of which 
increases with time and the other which decreases (Belton et al., 
2011a). In Section 5 we deduce the kinematics of the spin state 
and evaluate the magnitudes of M and E. In Section 6 we analyze 
spatially resolved images of the nucleus taken during the time 
interval Enc(ounter) - 1 h to Enc using position angle (PA) of the 
axis of minimum moment of inertia, i.e., the “long” axis, to deter- 
mine the direction of M and improved estimates of the intermedi- 
ate and short principal axes on the sky. In Section 7 we discuss the 
implications of the changing the estimated direction of the inter- 
mediate principal axis of inertia and of the changing spin state 
(and implied net torques) on conditions in the inner coma includ- 
ing the production of water vapor at the nucleus surface, the orien- 
tation of a plume of water vapor, and the time evolution of 
insolation at active areas. 


2. Relevant properties of the nucleus 

In Table 1 we collect various parameters that describe the ori- 
entation of the nucleus at the time of closest approach and are 
needed to specify the spin state. Most of these quantities have been 
noted earlier by Thomas et al. (2012), but we repeat them here for 
convenience. Also needed in this study is a competent model of the 
shape of the nucleus (Fig. 1 ). This model, again provided by Thomas 
et al. (2012), is available on a 2° longitude-latitude grid with its 


Table 1 

Relevant properties of the nucleus. The tabulated values are mostly from Thomas et al. (2012) and are collected here for convenience. Other values are from Thomas (private 
communication). The principal moments of inertia //, /„ and I s are related to the long (/), intermediate (i), and short (s) axes respectively. In this study we use the // and /, moment 
orientations to compute the direction of I s from their cross-product according to rule s = l xi (Samarasinha and A’Hearn, 1991 ). As a result, the direction of I s given in the table is 
the reverse of that given by Thomas et al. (2012), who use the convention s = ixl. The values for the moments of inertias are computed assuming that the interior mass 
distribution is homogeneous and are given per unit mass of the nucleus. Orbital information was obtained from the JPL HORIZONS web site where the spacecraft is referred to by 
“@-140”. 


Quantity 

Value 

Time of closest approach (UTC at comet) 

2010-11-04, 13:59:47.7 
JD 2455505.0831866 

Mean radius (m) 

580 ± 18 

Ratios of model principal moments of inertia 

him 2 ) 

5.302 ±0.4 x 10 4 

1,11s 

0.166 ±0.004 

I, ft 

0.979 ± 0.002 

Body moment orientations (latitude, E longitude) 

h 

89.73°, 207.56°E 

Ii 

0.05°, 106.14°E 

Is 

-0.25°, 196.14°E 

Coordinates of the center of mass (latitude, longitude, and radius vector (m)) 

+19.57°, 240.39°E, 5.0 

Direction of long axis (small end) at closest approach (RA, Dec; J2000) 

RA: 226.12 ±0.9° 
Dec: 39.37 ± 0.7° 

Direction of minimum moment at closest approach (RA, Dec; J2000) 

RA: 225.92 ± 0.9° 
Dec: 39.60 ±0.7° 

Direction of intermediate moment at closest approach (RA, Dec; J2000) [the value given is that provided by 

RA: 330.0° 

Thomas et al. (2012). In the solution for M found here, this direction had to changed to 302°, -16.5° (see 

Dec: 16.3° 

text)] 

(the on-sky uncertainty is ± 12° in a band normal to the 
direction towards the long axis) 

Surface area (m 2 ) 

5.24 x 10 6 

Volume (m 3 ) 

8.09 x 10 8 

Length of the long axis (m) 

2330 

Estimated mean density (kg/m 3 ) 

230-300 

Estimated mass of nucleus (kg) 

1 .9-2.4 x 10 11 
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Fig. 1 . Comparison of MRI images with model simulations of the nucleus of 103P at E-199.07s, E-0.67s and E+403.88s. The images are mvl0110413_5004000, 
mvlOl 10413 J5000002, and mvlOl 10414_5006046. The sub-solar region can be seen as a dark area in the “waist” region of the two left images. The shape of both the nucleus 
and terminator are simulated reasonably well even though, in this rendition, the scaling (not orientation) of the images and models were approximated by hand. The changing 
spatial resolution is evident in the images and is a reflection of the different ranges at which the images were taken. From left to right the ranges are: 2548, 694, and 5024 km. 


various radius vectors originating near the center of figure. We use 
this model, as shown in Fig. 1, after its transformation to a cylindri- 
cal coordinate system with the origin at the center of mass and 
coordinate positions referred to the principle axes of inertia. As ap- 
plied to this problem, we have transformed the model to a coarser 
grid mainly to reduce computing time and the sizes of the matrices 
involved (all the calculations in this paper were done in MathCad 
scripts [Mathsoft Engineering & Education, Inc.] or in Microsoft Ex- 
cel spreadsheets). The location of the center of mass assumes a 
homogeneous interior (see Section 3). 

Although our model is relatively coarse, it retains enough infor- 
mation on the shape to display its primary features and, in order to 
compute the time dependence of solar illumination at any location 
on the nucleus (Section 7), normals to the local surface in the 5°- 
grid model were computed at any point by using information from 
it and the surrounding four grid points. 

3. Overview of method and assumptions used to determine the 
spin state 

Our basic assumptions are that the mass distribution in the 
interior of the nucleus is homogeneous and the nucleus rotates 
as a rigid body. These assumptions were also made and justified 
for cometary nuclei in an earlier study of the spin of 9P/Tempel 1 
(Belton et al., 2011b), and we refer the reader to that work for de- 
tails. Ultimately, however, the justification rests on how well the 
model simulates the observed motion, and we will present an 
example in Section 6. With these assumptions, we can make use 
of the detailed theoretical description of torque-free asymmetric 
rigid-body spin developed in the Appendix of Samarasinha and 
A’Hearn (1991). We will not repeat the lengthy development of 
their analysis here except to note that we use precisely the same 
nomenclature, definitions, and system of Euler angles that are used 
in Samarasinha and A’Hearn’s work. When it is necessary to refer 
to one of their equations, we use the same numbering system as 
found in Samarasinha and A’Hearn followed by “(SA).” For exam- 
ple, the fundamental equation of motion ( dM/dt) ext = 0 is Eq. 
Al(SA). 

We begin with the estimation of the periodicities P^(t) and 
Pi//(t), which are the time averaged period of circulation of the long 
axis about M and the period of roll or oscillation about the long 


axis at time t respectively. These quantities are evaluated at differ- 
ent times during the encounter by dividing the comet’s 81.5d light- 
curve into sections. The subscripts </>, x// (as well as a third, 6), refer 
to the same Euler angles that are defined in Fig. Al of Samarasinha 
and A’Hearn (1991). </> and if refer to the rates at which the long 
axis circulates about M and the rate of roll or oscillation about 
the long axis, respectively. 6 is the angle subtended between the 
long axis and M, which oscillates with a period of P^/2 around a va- 
lue that we find to slowly decrease throughout the encounter. 

The fact that we find that and P^ vary secularly with time im- 
plies that the motion is not torque-free, and this finding requires 
that we make a further assumption, again to be justified a posteri- 
ori. We assume that the actual motion throughout a single spin cycle 
deviates at most by a small amount from torque free motion. With 
this assumption we can use the formalism of Samarasinha and 
A’Hearn (1991) to calculate the spin state at each instant of time 
with the relevant values of the evolving rotational momentum 
and energy. In this work the evolution of the spin state is, in effect, 
assumed to be a progression of torque-free states. 

Once P^(t) and P^(t) are known, we can determine which spin 
mode is appropriate and then estimate the magnitude of M(t) 
and the rotational energy, E(t) throughout the encounter. The Euler 
angles i f{t) and 0(t) can also be evaluated at this point. 

A first approximation to the space direction of M, which is as- 
sumed to be fixed in this study, is then estimated from i/^(0) and 
0(0), and the space directions of the long and intermediate princi- 
pal axes at closest approach from Thomas et al. (2012), where t = 0 
is the time of closest approach. This estimate, whose uncertainty is 
dominated by the considerable uncertainty in Thomas et al.’s 
(2012) estimate of the direction of the intermediate axis, is then 
improved by applying high-resolution observations of the chang- 
ing position angles (PA) of the nucleus, measured over as long a 
time interval as possible, as constraints during this time by mini- 
mizing variance between the model and the data. 

4. Determination of P^(t) and P^(t) from encounter light-curves 

Aperture photometry, centered on the nucleus, was performed 
on 32,049 calibrated images using the IDL “aper.pro” program with 
a 3-pixel radius aperture. The photometric reductions generally 
followed the procedures that were used in an earlier study of Co- 
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using the ANOVA technique (Schwarzenberg-Czerny, 1996) as 
implemented in the commercial Peranso astronomical software 
package (www.CBABelgium.com). We generally work in the fre- 
quency domain (cycles/d) rather than the time domain (periods), 
because it is much easier to perceive harmonic relationships be- 
tween individual spectral features in the power spectra in the fre- 
quency domain. The ANOVA software, which was used successfully 
in our study of 9P/Tempel 1, displays the power spectrum and 
computes the frequency of selected spectral features and their 
uncertainty. The method of estimating these uncertainties is that 
described in Schwarzenberg-Czerny (1991). In Fig. 3a and b we 
70 .$0 -so -40 -30 -20 -10 0 10 20 30 show the groups of data and the resulting power spectra. 

These spectra show considerable changes in the distribution of 
power amongst the spectral features (lines) as the encounter pro- 
gresses. In addition, there are systematic shifts in the peak periods 
of the lines. This is shown in Table 3 and in Fig. 4. There are two 
groups of features. One, Group 1, contains three lines, marked as 
Di, r> 3 , and v 4 in the figure, that generally decrease in frequency 
as the encounter progresses. The second group, Group 2, contains 
two lines (v 2 ,v 5 ) that generally increase in frequency as time pro- 
gresses. Evidently, these two groups are related to different phys- 
ical aspects of the rotation. Inspection of the numerical values of 
the frequency of the lines shows that the lines in each group are 
harmonically related, i.e., r>i ^ 2v 3 = 3d 4 and v 2 = 3r> 5 . This rela- 
70 -SO -30 -40 -30 -20 -10 0 10 20 30 tionship is illustrated as frequency differences that hover around 

Time from Encounter (d) zero in Table 4. 




Fig. 2. (Top) 103P light-curve using the MRI camera with the CLEAR1 filter. This 
data consist of 24,799 points and have been scaled to remove the inverse-square 
dependence on heliocentric distance and range. (Bottom) Same as the upper panel 
but with the background removed. The removal process was done in —10 day 
sections (see Table 2). 

met 9P/Tempel 1 (Belton et al., 2011b). A more complete descrip- 
tion of the overall reduction procedure in provided in Bodewits 
et al. (2012). Only MRI (Medium Resolution Instrument) images 
in the CLEAR1 filter, not compromised by cosmic ray hits, satura- 
tion, smear, or other evidence of corruption, were used. This left 
some 24,799 measurements to be used in the analysis (Table SI 
in Supplementary Online Data). The data cover 81.6 d during the 
period from Enc - 60 to Enc + 21.6 d. Two versions of the light- 
curve, one background subtracted, are shown in Fig. 2. The original 
data were adjusted to remove the dependence on heliocentric dis- 
tance and a presumed inverse square dependence on range. Since 
the signal is a mixture of resolved coma and unresolved nucleus, 
the latter step is not essential, but its inclusion does not adversely 
affect the analysis. The data were divided into seven groups on the 
basis of date rather than number of observations (Table 2). This 
grouping was done to better spread out the derived periodicities 
in time, even though the power spectra had very different signal- 
to-noise ratios. Each group was analyzed for latent periodicities 


Table 2 

Section of data used in the analysis of the light-curve. The median time is defined as 
the average time of all of the data in the block. The times are measured from the time 
of closest approach (Table 1). 


Data 

section 

Start 
time (d) 

Stop time 

(d) 

Median 
time (d) 

Number of 
points 

Sampling 

1 

-60.0209 

-48.6907 

-54.5876 

318 

Sparse 

2 

-48.6904 

-40.0217 

-44.3603 

684 

- 

3 

-33.9018 

-24.1916 

-29.0365 

3897 

- 

4 

-23.9018 

-14.2367 

-19.4743 

4061 

- 

5 

-13.9018 

-0.2067 

-7.8314 

4728 

- 

6 

0.1878 

11.9993 

7.0710 

10016 

Continuous 

7 

12.0003 

21.5698 

16.5403 

1095 

Continuous 


4A. Physical identification of the observed frequencies 

Delay-Doppler radar images (Harmon et al., 201 1 ) and a double- 
peaked nucleus light-curve (Meech et al., 2009, 2011) definitively 
relate observations of an increasing period of 16.6-18.2 h, i.e., a 
decreasing frequency from -1.45 to 1.32 c/d, with the circulation 
(rotation) of the long axis of an elongated nucleus about M in the 
interval between April-May, 2009 and during October, 2010. This 
rate is in good agreement with the magnitude of the change in 
Group 1 and provides a secure identification of the physical cause 
of t?i : the circulation of the long axis of this elongated body around 
M. r> 3 and v 4 are clearly sub-harmonics of Oi presumably stimu- 
lated by the locations of whatever source(s) of coma material are 
ultimately responsible the light-variations. 

The theory of rigid body rotation indicates that a second group 
of frequencies is related to either a roll or oscillation about the long 
axis (Samarasinha et al., 2004). Evidence for a second periodicity 
has already been seen indirectly in Earth-based observations. 
Microwave spectroscopy near encounter by Drahus et al. (2011) 
shows that the pattern of HCN production repeats every three 
—18.2 h cycles, and Knight and Schleicher (2011), studying coma 
structures, found that “the best rotational matches were found 
for images multiples of three cycles apart.” More recently, Waniak 
et al. (2012) have found that their images of the central coma also 
have a three-cycle repetition. In addition, a similar 3-cycle pattern 
can be seen in the post-encounter DIXI data, e.g., in the CLEAR1 fil- 
ter light-curve (Fig. 3b) for the 12 days immediately following clos- 
est approach. Evidently, the rotation of the nucleus brings the 
sources on the nucleus to essentially the same configuration with 
respect to the Sun every three-cycles, i.e., -every 54-55 h or with 
a frequency near 0.44 c/d. We can be reasonably assured of this re- 
sult if we identify v 2 , but not v 5 , as the frequency of roll or oscilla- 
tion about the long axis since, as seen in Table 5, [v^}l{v 2 ] averages 
1.5 throughout the encounter. v 5 is therefore a sub-harmonic of v 2 , 
i.e., v 2 = 3d 5 (Table 4). 

In summary, we identify Oi with the circulation of the long axis 
around M and v 2 with the roll or oscillation of the nucleus around 
its long axis. In order to make maximum use of the information in 
the power spectra and the precision of the harmonic relationships 
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Fig. 3. (a) (Left side) Sections 1 through 4 of the light-curve with Section 1 at the top. (Right side) Periodograms (power spectra) computed using the ANOVA method showing 
the evolution of the spectral content, (b) As in 3(a) - but for data Sections 5-7. 


seen in Table 4, we combine the information of all the lines within 
each group to form a weighted average {bi }(t) and { v 2 }{t ) at each 
time step. Thus, P^ = and P ^ = {l 2 } _1 . The resulting values 

are collected in Table 5, and the time-evolutions of {ui}(t), { v 2 }{t ), 
P^(t) and P^(t) around the time of encounter are shown in Fig. 5. 
We emphasize that the relationship between {bi} and {b 2 } is 
approximate and varies slowly throughout the encounter. It is 
not, apparently, a constant of the motion or indicative of a reso- 


nance. The period of circulation of the long axis is seen to steadily 
increase through perihelion passage in an approximately linear 
fashion at a rate of ~1.3 ± 0.2 min/d in agreement with what was 
found in the Earth-based observations of Drahus et al. (2011) and 
Knight and Schleicher (2011). The evolution of the roll around 
the long axis is more complex, with the angular rate apparently 
slowing down at first until about 30 d before perihelion and then 
speeding up through perihelion passage. However, the reality of 
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Fig. 3 ( continued ) 


this early slow down is questionable since the evidence rests on 
the value of a single point. The rate of change of roll period is 
approximately -4.4 min/d at perihelion. There is a hint in the data 
that the rate of change of roll is beginning to stabilize some 20 d 
after perihelion, but this conclusion is obviously very uncertain. 


principal moments of inertia //, I it I s (Table 1). Samarasinha and 
A’Hearn (1991, Appendix) have developed two inequalities, where 
the value of the ratio P^/P^ relative to functions of the principal 
moments of inertias, determines whether the spin mode is a 
LAM or SAM. 

For a LAM we must have: 


5. Kinematics of 103 P/Hartley 2*s spin state 

We look at the kinematics of 103P’s spin in two parts. In the first 
we determine the mode in which the nucleus is spinning and com- 
pute M (the magnitude of M), E and other parameters of the spin 
state based on the above values of P^ and P^. In the second we 
determine our best estimation of the space orientation of the spin 
state using PAs of the long axis of the nucleus measured in pre- 
encounter observations. 



IJs 

( li-h)(Is-h ) 


-1 


(A52(SA)) 


and, for a SAM: 




IJi 

m 



(A78(SA)) 


( 1 ) 


( 2 ) 


5.1. The mode of spin 

The question of what spin mode the nucleus finds itself in dur- 
ing the encounter depends not only on P^ and P^ but also on the 


In Fig. 6 we apply these inequalities to the situation for 103P. 
Throughout the encounter the spin state is seen to lie firmly in 
the LAM region, with the decrease in this ratio with time indicating 
an increase in the level of excitation as the encounter progresses. 
Also shown in the figure are the effects of the level of uncertainty 
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Table 3 

Evolution of line periodicities (h) through the encounter. 


Time from encounter (d) 

Periodicies (h) 





ib 

02 

v 3 

v 4 

05 

-54.5877 

17.05 ±0.23 

27.42 ± 0.37 

34.28 ± 0.57 

51.57 ±0.86 

87.09 ±2.59 

-44.3603 

17.10 ±0.78 

28.28 ±0.80 

34.28 ±1.51 

53.69 ±5.76 

86.01 ± 6.07 

-29.0365 

17.83 ±0.59 

28.62 ±0.39 

35.61 ±0.95 

58.16 ±2.24 

86.14 ±1.77 

-19.4743 

18.02 ±0.06 

28.28 ±0.13 

36.11 ±0.16 

54.46 ±0.14 

84.66 ± 0.32 

-7.8315 

18.12 ±0.05 

27.13 ±0.28 

36.17 ±0.08 

54.36 ± 0.07 

76.29 ±6.16 

7.0711 

18.63 ±0.02 

26.32 ±0.26 

37.26 ± 0.03 

55.76 ±0.04 

79.00 ± 0.68 

16.5404 

18.61 ±0.09 

26.13 ±0.23 

37.28 ±0.24 

55.92 ±0.30 

78.19 ±2.09 
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Fig. 4. (Left) Evolution of the frequencies of the five most prominent lines in the power spectra of the seven blocks of data. The lines are labeled in order. r>i through v 5 . 
However, for clarity of presentation, we have labeled only the first and last lines in the Figure. (Right) Evolution of the periods of the five lines. In both parts of the figure, two 
groups of lines can be identified: (n lf v 3 , v 4 ) and (v 2 , o 5 ). In the first group the frequencies of the lines decrease with time through perihelion passage. In the second the 
frequencies increase with time. Perihelion passage (JD 2455497.756967) is marked as a vertical dotted line. 


Table 4 

Frequency (c/d) differences between base line frequency and their sub-harmonics in the two groups of lines. In most cases the magnitude of the differences are substantially 
smaller than the uncertainty in the base frequency. 


Group 1 




Group 2 



lb 

Uncertainty 

2o 3 - r>i 

3l> 4 - r>i 

02 

Uncertainty 

3d 5 - v 2 

1.28932 

0.00592 

-0.0016 

-0.00178 

0.91844 

0.00824 

0.00235 

1.28846 

0.00114 

—IE— 04 

0.00283 

0.91179 

0.00902 

-0.00042 

1.32436 

0.00344 

0.00264 

8E-05 

0.88457 

0.00924 

0.05923 

1.33204 

0.00464 

-0.00272 

-0.01012 

0.84868 

0.00396 

0.00182 

1.34592 

0.0448 

0.00196 

-0.10794 

0.83864 

0.01142 

-0.00281 

1.40384 

0.0642 

-0.00356 

-0.06287 

0.84852 

0.02396 

-0.0114 

1.40737 

0.01917 

-0.00697 

-0.01123 

0.87541 

0.01194 

-0.04867 


in the moments of inertia. Only if there were significant large-scale 
axial inhomogeneities in the internal mass distribution in the nu- 
cleus would it be possible for the nucleus to be in a SAM state, 
i.e., even though the nucleus is close to axial symmetry, the mo- 
ments Ii and I s would have to have radically different values. Further 
investigation on this point is beyond the scope of this work, and we 
conclude that the nucleus of 103P spins in the long-axis mode with 
the level of excitation increasing through the time of closest 
approach. 


5.2. Determination of rotational angular momentum and energy of the 
spin 


As shown in the Appendix to Samarasinha and A’Hearn (1991) 
we can set up definitions for the observed P ^ and P^ in terms of 
M and E, i.e., 


(P-,)=2n. 


r p f / 2 

(2/P*) / < 

Jo 


ij)dt 


(SA not numbered) 


(3) 
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Table 5 

Weighted estimates of the period (frequency) of circulation of the long axis about M and the roll period (frequency) throughout the encounter. 


Time from encounter (d) 

{»i} (c/d) 

Aui (±) 

P* (h) 

AP<j) (±) 

f» 2 } (c/d) 

Al> 2 (±) 

P*( h) 

AP\|/ (±) 

{vMv 2 } 

-54.58 

1.4022 

0.0125 

17.12 

0.15 

0.8661 

0.0107 

27.71 

0.34 

1.618 

-44.36 

1.3967 

0.0425 

17.18 

0.52 

0.8469 

0.0222 

28.34 

0.74 

1.649 

-29.03 

1.3193 

0.0242 

18.19 

0.33 

0.8378 

0.0095 

28.65 

0.33 

1.574 

-19.47 

1.3260 

0.0025 

18.10 

0.03 

0.8498 

0.0025 

28.24 

0.08 

1.560 

-7.83 

1.3249 

0.0013 

18.11 

0.02 

0.8854 

0.0092 

27.11 

0.28 

1.496 

7.07 

1.2897 

0.0006 

18.61 

0.01 

0.9116 

0.0059 

26.33 

0.17 

1.414 

16.54 

1.2884 

0.0039 

18.63 

0.06 

0.9187 

0.0078 

26.12 

0.22 

1.402 




Time from encounter (d) 

Fig. 5. Evolution of weighted estimates of spin frequencies (Right) and periods (Left). The filled circles represent the circulation of the long axis around M, and the filled 
triangles the roll around the long axis. The vertical dotted line marks the time of perihelion passage. The reality of the sinuous behavior of the roll, which appears to slow 
down at first and then speed up through perihelion, is uncertain, as it depends on a single point. Only the trend through perihelion is secure. 


where 


</> = M 


(Ii-h) + (I s -Ii)*sn 2 z)' 
/,(/| - /,) + /,(/, - /|)OT2T 


(A38(SA)) 


( 4 ) 


and 


Pit = 4 


w. 


^2 


n/2 


du 


1,2 sin 2 u 


(A45(SA)) (5) 


where 


k 2 = 


and 




(A32(SA)) 


( 6 ) 


T = t* 


2 E(J, 


W. 


(A31(SA)) 


( 7 ) 


In Eq. (3) the integration is taken over P^/2, because it is one of the 
peculiarities of rigid body motion that is periodic over this inter- 
val. Eqs. (3) and (5) are non-linear but can be solved numerically if 
we have good first approximations for M. and E. We get a first 
approximation to the magnitude of M, i.e., M, from 

(P») mi„-^ (A49(SA)) (8) 

where the arrow indicates the limiting value, and 
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Fig. 6. Determination of the mode of spin. The ratio falls entirely in the area 
where LAM states are allowed (above the dash-dot line defined by Eq. (1)). The 
(unbroken) horizontal line, shown with its uncertainty level, is defined by Eq. (2). 
Only above this line are SAM states allowed in this system. 
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(Pi) 


max 


2nIiVUh 

M[VIJ~ 5 -y/(I^Ii)(Is-Il)] 


(A50(SA)) 


( 9 ) 


by putting the observed value of equal to the functional form of 
((^)min + (^)max)/ 2 an d solving for M. Using this value and the 
observed value for P ^ in Eq. (5), we compute a first approximation 
to E. With these estimates as starting values we can then solve Eqs. 
(3) and (5) iteratively for precise values of E and M. In this study 
we use the MathCad “solve block” Find function to perform this 
iteration. The results of this process are collected in Table 6, which 
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Table 6 

Values of M, E, and k 2 for 103P and their variations throughout the encounter, k 2 (Eq (6)) is an inverse measure of the level of excitation in the spin state. As the encounter 
progresses through perihelion, the momentum and energy both decrease, however the level of excitation increases. Columns 3, 5, 7 are estimates of the uncertainties in the 
various quantities. 


Time from encounter (d) 

M (m 2 /s) 

AM (±) 

E (m 2 /s 2 ) 

A E(±) 

k 2 

A (k 2 ) (±) 

-44.36 

3.208E+01 

6.90E-01 

1.750E-03 

7.03E-05 

2.454E-01 

0.201 

-29.04 

3.030E+01 

3.92E-01 

1.572E— 03 

3.77E-05 

2.262E— 01 

0.078 

-19.47 

3.046E+01 

3.98E-02 

1.590E— 03 

3.87E-06 

2.225E-01 

0.006 

-7.83 

3.043E+01 

2.13E-02 

1.598E-03 

2.82E-06 

2.063E— 01 

0.003 

0 

2.996E+01 

1.48E-01 

1.557E— 03 

1.40E-05 

1.954E— 01 

0.020 

7.071 

2.963E+01 

9.55E-03 

1.529E-03 

1.57E-06 

1.864E-01 

0.001 

16.54 

2.959E+01 

6.39E-02 

1.529E— 03 

6.23E-06 

1.833E— 01 

0.007 


again shows that the excitation of the spin state is actually increas- 
ing through perihelion even though the total rotational momen- 
tum and energy are steadily decreasing. The trends in M, E , and 
k 2 (Eq. (6), an inverse measure of excitation) are shown in Figs. 7 
and 8. 

5.3. Determination of other spin parameters at closest approach 

Once M and E are determined other parameters that are a func- 
tion of time can be calculated according to Eqs. A31(SA)-A38(SA). 
These parameters include two of the Euler angles (i J/ and 0), the 
components of the total spin vector, and the (variable) rate at 
which the long axis circulates around M. We leave the determina- 
tion of most of these quantities to the following section where, we 
estimate the orientation of M and the nucleus in space. However, 
there are aspects of the motion that depend on the degree of exci- 
tation and exhibit secular change throughout the encounter, which 
we evaluate here. These quantities include the maximum and min- 
imum values of 9 , the Euler angle between the long axis and M, and 
the maximum and minimum values of the roll rate, \j/, and the cir- 
culation rate, </>. The secular evolution of these quantities through- 
out the encounter is shown in Table 7. The mean value of 6 steadily 
decreases from 82.8° to 81.7° (another indication of increasing 
excitation), while the full amplitude of the associated “nodding” 
motion, which has a period of P^/2, decreases from ~1° to 0.8°. 
The rate of roll and circulation are shown to vary within each spin 
cycle by 14-10% and by about 2% respectively. 


6. Determination of the orientation of the spin state in space 

As we noted in the introduction, there is little consensus in the 
various studies of coma morphology on the space direction of M for 
103P. This scatter is shown in Table 8. In this work we have the 
benefit of the radar observations of Harmon et al. (2011), which 
show definitively that the circulation of the long axis is in the pro- 
grade sense, and the determination of the space directions of the 
minimum and intermediate axes at the time of closest encounter 
(Thomas et al., 2012). We note at this point that the sense of roll 
around the long axis (anti-clockwise looking at the small-end of 
the nucleus along the long axis) is fixed once the sense of the mo- 
tion about M is determined (Ames and Murnaghan, 1929). To- 
gether with the Euler angles, specifically 9(0) and it is a 
matter of simple spherical trigonometry to estimate the space 
direction of M at that time. 

The Euler angles at time t are calculated from: 


m = atan2 j , k), yC-^cnr(t, k) j (A36(SA)) 

( 10 ) 

where the order of the parameters in the atan2 function is as de- 
fined in Fortran, 



Fig. 7. (Top) The evolution of M through the encounter. The rotational momentum 
is seen to decrease steadily through perihelion passage. (Bottom) The rotational 
energy also decreases through perihelion passage, which is marked with a vertical 
dotted line. 
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Fig. 8. The index k 2 (Eq. (6), an inverse measure of the degree of excitation) shows 
that the spin state moves to a higher excitation level as the encounter passes 
through perihelion, which is marked with a vertical dotted line. 


604 


M.J.S. Belton et all Icarus 222 (2013) 595-609 


Table 7 

Secular variation of 9, and </>, x Jy motions throughout the encounter. 


Time from encounter (d) 

9 min (°) 

9max (°) 

•Amin (C/d) 

•Amax (c/d) 

•Amin (C/d) 

•Amax (c/d) 

-44.36 

82.47 

83.46 

0.7866 

0.9095 

1.3812 

1.4109 

-29.04 

82.16 

83.11 

0.7833 

0.8943 

1.3047 

1.3327 

-19.47 

82.10 

83.04 

0.7955 

0.9060 

1.3114 

1.3396 

-7.83 

81.80 

82.70 

0.8333 

0.9394 

1.3104 

1.3385 

0 

81.57 

82.45 

0.8484 

0.9498 

1.2899 

1.3176 

7.07 

81.37 

82.23 

0.8633 

0.9612 

1.2756 

1.3030 

16.54 

81.30 

82.15 

0.8710 

0.9679 

1.2743 

1.3016 


0(t) = COS 




I,(I s -M 2 /2E) 


(M Z /2E)(I S — If) 


(A37(SA)) 


( 11 ) 


and 


= f 4>dt. (12) 

Jo 

0(t) is given in Eq. (4), and snr, cm and dm are Jacobian elliptic 
functions. 0(0) is determined once the direction of M has been cal- 
culated from \j/( 0), 0(0) and the directions of the minimum and 
intermediate axes of moment of inertia. In order to calculate 0(0) 
we also choose to orient the X-axis of the angular momentum coor- 
dinate system (Fig. A1 in Samarasinha and A’Hearn (1991)) into the 
meridian plane of the Equatorial coordinate system that contains M. 
At closest approach we find 0(0) = 241.6°, \j/( 0) = 0° and 
0(0) = 81.60°, and given the direction of the inertia axes in Table 1, 
we find our first approximation to the space direction of M as (RA, 
Dec; J2000) = 83.18°, +53.39°. 

The uncertainty in this approximation to the direction of M is 
substantial, since it originates in the combined uncertainties in 
our estimates of the moments of inertia, M, E, and the directions 
of the long and intermediate axes. The uncertainty in the latter is 
particularly severe, since according to Thomas et al. (2012) it 
amounts to ±12° in a band normal to the direction of f (the long 
and intermediate axes are, by definition, orthogonal), so the most 
likely values of the RA and Dec of M lie on a highly elongated area 
across the sky centered on the first approximation as shown in 
Fig. 9. Deconvolved DIXI high-resolution images (Lindler et al., 
2012) provide a powerful constraint here by providing, through 
measurements the PAs of the long axis over an arc of ~20° on 
the sky during ~1 h of the time before closest approach, a measure 
of the orientation of the track of the rotating long axis on the sky. 

In Fig. 10 we show a sample of deconvolved HRI (High-Resolu- 
tion Imager) images taken in the period Enc - 1 h to Enc + 1 h. In 
our determination of the PA of the long axis we found that only 
images in the interval Enc - 1 h to Enc were found suitable owing 
to the peculiar illumination of the terminator region in the post- 
encounter images (Fig. 10). In Table 9 we list the images that were 
used and our estimated PAs. These estimates are the mean of three 
independent sets of measurements by two individuals. Based on 
the dispersion in these measurements, we assign an uncertainty 
of ±1.7° to the PAs. The results are plotted in Fig. 11, which show 
how the projected orientation of the nucleus changes on approach, 


but still well in advance of the time of closest approach when the 
changes are dominated by the spacecraft orbit. Also shown in 
Fig. 1 1 is a curve that shows the predicted change in PA based on 
the first approximation of the direction of M (X = 0°). A clear diver- 
gence exists between the predictions of this first approximation 
and the observed PAs. The origin of this discrepancy is, as shown 
in the figure, the large uncertainty associated with the direction 
of the intermediate principal axis of inertia. In the figure we test 
the result of moving the direction of the intermediate axis along 
its arc of uncertainty by X° (Fig. 9). An excellent account of the ob- 
served PAs is found for X = 43° (Fig. 11), where M is directed to- 
wards (RA, Dec; J2000) = 8±4°, +54 ±1°. This our best estimate 
of the direction of M (equivalent to an obliquity of 48 ± 1°), which 
determines the orientation of the state. Fig. 12 shows change of 
variance with X in this determination. 

The space direction of the instantaneous spin vector at this time 
is (RA, Dec; J2000) = 300 ± 3°, +67 ± 1.3°. The “adjusted” direction 
of the intermediate principal axis of inertia is (RA, Dec; 
J2000) = 301.8°, -16.5°, some 33° away (on the sky) from the Tho- 
mas et al. (2012) estimate, and the corresponding adjustments to 
the initial Euler angles gives 0(0) = 298.76°, = 0°, and 

0(0) = 81.60° . The direction of the short principal axis of inertia 
must therefore be moved to 14.1°, +45.8° for consistency. The 
instantaneous spin vector is inclined at an angle of 33.2 ± 1.3° to 
M, and the instantaneous spin period is 14.1 ± 0.3 h at closest ap- 
proach. Fig. 13 shows the alignment of the model and observation 
for the last resolved HRI image after encounter for our best esti- 
mate of the spin state. 


7. Discussion and conclusions 

7.1. Summary of the derived spin state and its evolution 

At the time of closest approach the nucleus is found to spin in 
the long-axis mode (LAM) with the long axis circulating around 
M with a period near 18.40 ± 0.13 h, tilted to M near an angle of 
81.2 ± 0.6°, while simultaneously rolling around the long axis with 
a period of 26.72 ± 0.06 h. The motion is in the prograde sense. As 
an alternative description, the instantaneous spin vector moves 
through the body of the nucleus circulating around M every 
1 8.40 ± 0.1 3 h, while tilted near an angle of 33.2 ± 1 .3° and rotating 
with an instantaneous period near 14.1 ± 0.3 h. These periodicities 
and tilts are not constant during the motion and vary by a few 


Table 8 

Estimates of the direction of M. 


Reference 

Observation type 

RA (°) 

Dec (°) 

Notes 

Knight and Schleicher (2011) 

Coma features 

257 

+67 

- 

Samarasinha et al. (2010) 

Coma features 

345 ± 20 

-15 ± 20 

- 

Waniak et al. (2012) 

Coma features 

122 

+16 

- 

Harmon et al. (2011) 

Radar 

17 

47 

Consistency 

P.A Taylor (Private Communication) 

Radar 

347 ± 20 

+27 ± 20 

Preliminary 

Belton et al. (this work) 

Nucleus imaging 

8 ± 4 

54 ±1 

- 
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Fig. 9. Estimates of the direction of M. The direction found in this work is marked 
with an open rectangle. The first approximation to this direction, based on the 
original Thomas et al. (2012) direction of the intermediate principal axis of inertia, 
is marked with a filled circle. The line joining them on the sky is the locus of 
possible directions given the uncertainty in the latter (see text). The direction found 
by Knight and Schleicher (2011) is marked with an asterisk; that by Samarasinha 
et al. (2010) by an open diamond; that by Waniak et al. (2012) by an open circle; 
and a radar based direction communicated by P.A. Taylor by an open triangle. 
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Fig. 10. Sample of three deconvolved pre-encounter HRI images (Left column) and 
three post-encounter images (Right column) that cover the interval Enc - 1 h to 
Enc + 1 h. The changing position angles (PA) of the long axis, which easily be seen, in 
1 1 such pre-encounter images were used to determine the direction of M. Not all of 
the change in PA seen in this illustration is due to rotation. A fraction is due to the 
motion of the spacecraft relative to the comet. The post-encounter images, which 
have a more complex terminator structure, were not measured for PA but were 
used to check the result (Fig. 13). 

percent during a complete cycle. The direction of M on the sky (as- 
sumed constant in this study) was found to be (RA, Dec; 
J2000) = 8 ± 4°, +54 ± 1°; however, the direction of the intermedi- 
ate principal axis of inertia at closest approach estimated by Tho- 
mas et al. (2012) had to be adjusted by 43° along its arc of 
uncertainty to (RA, Dec; J2000) = 301.8°, -16.5° to satisfy observa- 
tions of the projected long axis during the encounter. 

The spin state was found to evolve during the encounter with 
the magnitudes of M and E slowly decreasing through perihelion 
passage at rates of -4.4 x 10 -7 m 2 /s 2 and 4.0 x 10 -11 m 2 /s 3 (per 
unit mass of the nucleus), respectively. However, even as the total 
energy decreased, the degree of excitation in the motion increased 
(Fig. 6). The average angle between the long axis and M decreased 
and as the average circulation period lengthened the average roll 
period decreased. To better illustrate the variability of the spin 
state we show in Fig. 14 predictions of the changing space direc- 
tion of the instantaneous spin vector and, in Fig. 15, the changing 
instantaneous spin period. 


The large displacement of the direction of the intermediate axis 
from the estimate by Thomas is perhaps not surprising given the 
near axial symmetry of the nucleus (in true axial symmetry there 
is no preferred direction for this axis) and possible imperfections 
in the shape model. It may indicate that some level of inhomgene- 
ity exists in the nucleus interior, but this possibility is hard to 
quantify other than to say that it would be small. 

In their analysis of the evolution of comet nucleus rotation, 
Neishtadt et al. (2002) find that axially symmetric comet nuclei 
with only a few active areas should tend to gain angular momen- 
tum, and align M with the direction toward perihelion. In the case 
of 103P, a comet nucleus with similar properties, M is found to be 
some 46° from the direction of perihelion, It is also losing rota- 
tional angular momentum. 

7.2. Relationship to other observations 

The somewhat complicated spin state described above clearly 
requires confirmation from observations other than those from 
which it was derived. These confirmations come from independent 
Earth-based observations, and we have already noted that the 
three 18.4 h cycles seen in several ground-based measurements 
are already built-into our model with the choice of 
P^( 0) = 26.72 h. The long axis of the nucleus circulates three times 
for every two roll periods to bring the nucleus back to essentially 
the same aspect with respect to the Sun. 

The secular evolution of the circulation period found in our spin 
state throughout the encounter is also found in ground-based 
observations. This confirmation is shown in Fig. 16 where the pres- 
ent results embed well with the ground-based record. 

It is in the direction of M where there is a considerable diver- 
gence of opinion among ground-based studies. However, we note 
that Knight and Schleicher (201 1 ) found that “ . . .one of the stron- 
gest constraints on the nominal rotation axis (i.e., the total angular 
momentum vector) is the side-on, corkscrew-like appearance of 
both the northern and southern features in early November. Not 
only do these features provide the best determined position angle 
through-out the apparition, yielding a great circle of deprojected 
pole solutions, but there is also the least amount of overlap be- 
tween the two features, directly implying that we are observing 
from close to the comet’s equator.” In their Fig. 1, the image taken 
at UT 2010 Nov. 11 09:39 shows this configuration very clearly. 
The angle between our direction for M and the line of sight from 
the Earth at this time is 94.5 ± 2.3°, i.e., the two directions are 
essentially orthogonal and so provide strong support for our result. 
The radar results of Harmon et al. (2011) also provide some sup- 
port, in that they find that the delay-Doppler images can be well 
explained by an unexcited rotation about a pole of RA=17°, 
Dec = 47°. Our spin state is but modestly excited (0~ 81° rather 
that 90°), and their assumed pole is only 9.0 ± 1.5° from our result. 
Similarly Harmon’s observations appear to be consistent with a 
simple rotation state with a spin axis in the region RA, 
Dec = (347, +27) ± 20° (P.A. Taylor, Personal communication), 
which partially overlaps our result. We conclude that Earth-based 
imaging and radar are consistent with our derived pole and there- 
fore provide a considerable degree of confidence in it. 

7.3. Torques 

A basic assumption in this study is that the direction of M can 
be assumed fixed during the encounter. It is obviously essential 
to check this assumption a posteriori using the measured secular 
rate of change in the magnitude of M through the encounter, i.e., 
using dM/dt= -4.39 x 10 _7 m 2 /s 2 per unit mass, which allows a 
measure of the net torque acting on the nucleus. With a total mass 
of 1. 9-2.4 x 10 11 kg (Table 1), we find a net torque of 
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Table 9 

Observations of the position angles (PA) of the long axis and other ancillary information in deconvolved HRI pre-encounter images used to determine the direction of M. In 
column 3, N. clock is the angle between the top of the camera frame and the projection of the North direction measured clockwise. 


Image 

Mid-exp. JD-2455505 

N. clock (°) 

PA (°) 

Range (km) 

Phase (°) 

hv0342 1 45836_5000005 _001 _rr 

0.0430539 

73.65 

44.75 

42731.44 

85.9 

hv0342 1 46502_5000022_001_rr 

0.0507665 

73.64 

49.46 

34522.65 

85.9 

hv0342 1 47067_5000065_001_rr 

0.0573033 

73.63 

52.08 

27565.93 

85.9 

hv0342147487_5000095_001_rr 

0.0621645 
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Time before encounter (min) 

Fig. 11. Plot of the PA’s (filled circles) measured on 11 pre-encounter deconvolved 
HRI images as a function of time from encounter. The curves represent the 
predicted run of PA based on directions of M that correspond to assumed directions 
of the intermediate principal axis of inertia that differ by X° from that found by 
Thomas et al. (2012). These directions fall on an arc of uncertainty that is normal to 
the direction of the long axis (see text). The coordinates assumed for M fall on the 
arc (line) shown in Fig. 9. The line marked X = 0°, which is well displaced from the 
observed points corresponds to the intermediate principal axis as found by Thomas 
et al. (2012). The curve for X = 43° provides the best accounting of the observations 
as indicated in Fig. 12. For this curve the direction of the intermediate principal axis 
of inertia is RA, Dec = 302°, -16.5°. 



Fig. 12. Dependence of variance (between model and observed PAs) on X. The 
minimum falls at X = 43°. 


~9.5 x 10 4 kgm 2 /s 2 was operating near perihelion. We can get a 
rough estimate of the upper limit to the precession rate for M by 
imagining that a torque of this magnitude is entirely applied 
orthogonally to M. Under these circumstances the precession rate 



Fig. 13. Comparison of the predicted model orientation and terminator with the 
last post-encounter deconvolved HRI image (hv0342150580-5007045_001_rr; mid- 
exposure at JD 2455505.09796810 or E+21.29 min) for our best estimate of the spin 
state. The image size has been adjusted but not the orientation in the image frame. 



Right Ascension (deg) 


Fig. 14. Sky tracks of the long axis (dotted line) and instantaneous spin vector 
(continuous line) over 100 h centered on the time of closest approach. The cross 
denotes our best estimate of the direction of M. Note that the spin vector does not 
follow a closed curve on the sky. 

is ~Torque/M = 1.5 x 10 -8 radians/s suggesting a maximum dis- 
placement of ~6° throughout the 81.6 d of the encounter. This 
upper limit is close to the estimated uncertainty in our determina- 
tion of the direction of M giving some confidence in our original 
assumption that the direction of M can be considered fixed during 
the encounter. 
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Fig. 15. The variability of the instantaneous spin rate over 100 h centered on the 
time of closest approach. 



Fig. 16. Observed values of (the period of circulation of the long axis around M) 
as a function of time. DIXI results are shown as filled circles. Other observations are: 
open circles - Knight and Schleicher (2011); open triangles - Samarasinha et al. 
(2011); open square - Meech et al. (2009); asterisk - Harmon et al. (2011); filled 
triangle - Drahus et al. (2011): cross - Waniak et al. (2012). Perihelion is marked 
with a vertical dotted line. The trendline is fitted, unweighted, to all of the 
observations except those of Meech et al. (2009). 



Fig. 17. High-resolution (8.5 m/pixel) MRI image (mv0342149264_5004032) was 
taken at JD 2455505.0827290 and shows where the evening and morning 
terminators meet on the big end of the nucleus. Images prior to this time are 
dominated by the morning terminator and those following by the evening 
terminator. In the image the evening part of the terminator is to the left. 


Belton (2012) has already discussed the possible relationship 
between the net torque and the gas (essentially water vapor; C0 2 
has small torque potential because it is a minor constituent and 
its direction of outflow is roughly parallel to the long axis) produc- 
tion from the nucleus using a preliminary estimate of the net tor- 
que. Since that time, we have increased our estimate of the torque 
by 27%. We will not repeat the argument in Belton (2012) except to 
note that with the improved estimate the surface of 103P evidently 
only contributes 0.8-6% of the total water vapor loss directly. The 
source of the bulk of the water lost by the comet apparently comes 
from icy particles in the coma that, to the first order, do not con- 
tribute to the net torque on the nucleus. 

7.4. Some applications 

7.4 A. Determination of the morning and evening terminators 

Given the shortness of the encounter relative to the spin period- 
icities, it is not obvious from the images or the shape model alone 
whether the terminator that is visible in an image is the morning 
or evening terminator. This distinction is determined by the nature 
of the spin. We have used our model of the spin state to resolve this 
problem. Some experimentation shows that the terminator visible 
on approach is the morning terminator and that post-encounter 
exhibits the evening terminator. In images at closest approach, it 
is the evening terminator that is visible. We find that the transition 
takes place near JD 2455505.08278; images before this time show 
the morning terminator and those after show the evening termina- 
tor over most of the nucleus. Close to the time of transition both 
terminators can be seen simultaneously at the big end of the nu- 
cleus. The high-resolution MRI image (mv5004032) shown in 
Fig. 1 of A’Hearn et al. (2011), and reproduced here as Fig. 17, 
was taken at JD 2455505.0827290 and shows both terminators 
simultaneously in the region of the big end of the nucleus. 

7.4.2. The evolution of insolation at specific regions of the nucleus 
We have computed the variation of cosine of the angle of the 
Sun direction with the normal (called here the “insolation”) over 
specific points of the nucleus to clarify our understanding of the 
time history of this driver for activity. In Figs. 18-20, we show 
the variability of the insolation with time over a 20 day period cen- 
tered on the time of closest approach over the tip of the small end, 
a point near the origin of one of the C0 2 jets at the small end (20°E, 
+79°), and a point near active areas seen in the evening terminator 
(120°E, -23°) respectively. At the tip of the small end the 18.4 h 
periodicity dominates as expected. The effect of the roll period is 
dramatic even for the high-latitude off-axis jets and certainly for 
the near equatorial region. The 3-cycle periodicity can be clearly 
seen in the latter two curves. 



Time from closest approach (days) 

Fig. 18. Time variability of the insolation (defined as the cosine of the angle 
between the normal to the surface and the direction to the Sun) over the tip of the 
small end of the nucleus from -10 to 10 days of encounter. The 18.4 h periodicity is 
easily seen. 



608 


M.J.S. Belton et all Icarus 222 (2013) 595-609 


Insolation factor for (Lat, Long) = (+79, 20E) 



Fig. 19. Time variability of the insolation over (Lat., Long.) = +79°, 20°E a point in 
the middle of the cluster of C0 2 jets seen near the small end of the nucleus -10 to 
10 days of encounter. The insolation is strongly modulated by the roll period and 
gives rise to the 3-cycle periodicity that can easily be seen. 


aT Insolation factor for (Lat. Long) = (-23, 120E) 

c ■ 
iS 

5 0.8 

"D 

£ 0.6- 
n 0.4- 
o 0.2 

in 


if 


1 

L 


1 





1 

1 


I 

1 




o 0 
O -10 

i -8 -6 -4 -2 

0 2 4 

6 

8 11 


Time from closest approach (d) 


Fig. 20. Time variability of the insolation over (-23°, 120°E) a region near the 
strong jets seen on the evening terminator from -10 to 10 days of encounter. 


7.4.3. The source of water vapor plume 

One of the remarkable results of the DIXI mission (A’Hearn et al. 
(201 1 ); see also Feaga et al., 201 1 ) was the discovery of a plume of 
water vapor emanating from the region of the waist of the nucleus. 
Fig. 21 shows the orientation of the nucleus and terminator in a di- 
rect comparison with A’Hearn et al.’s (2011) original rendering of 
the water vapor data. At the time these observations were made 
our model predicts that there were two sub-solar regions (defined 
as regions with solar zenith angle <24°) on the nucleus (because of 
its peculiar shape), one at the waist slightly toward the big end, 
and the other near the tip of the small end of the nucleus. We have 
included in the figure eight normals to the surface, four each from 
the boundaries that define each region and, as can be seen, one of 
the sub-solar regions, that near the waist, appears to be related to 
the source of the water vapor; however, the plume does not seem 
to propagate along the direction of the normals to the surface as 
might be expected if the plume issued as a jet from a vent in the 
surface. The normals to the surface are offset by ~27°. The direc- 
tion of the plume appears to be closer to a line that is orthogonal 
to the long axis and that passes through the sub-solar point (dotted 
line in Fig. 21), however, even this direction still diverges by ~15°. 
A possible explanation for this is that the direction of the plume is 
under the control of the local gravitational field, since Thomas et al. 
(2012) have shown that gravitational potential surfaces along a 
meridian flattens out above the waist region for bulk densities that 
are in the expected range for comets. The plume may simply be 
warm water vapor released in relative abundance in the sub-solar 
region that is buoyant in the relatively dense, i.e., collisional, ambi- 
ent atmosphere expected in the innermost coma and simply con- 
verts upwards. 

The other peculiarity that we would note is the lack of water va- 
por emanating from the sub-solar region at the small end 
(although, according to A’Hearn et al. (201 1 ), there is copious emis- 



Fig. 21. Comparison of the water plume found emanating from the waist of 103P/ 
Hartley 2 by A’Hearn et al. (201 1 ) with the predicted orientation of the nucleus and 
terminator. The prediction is for JD 2455505.0882922 in the middle of the IR scan at 
E+7 min. Also shown are normals from points of the surface that bound the two 
sub-solar regions (defined as regions with solar zenith angle <24°) that existed on 
the surface at the time. The base the water vapor plume is consistent with its source 
being the sub-solar region near the waist. However, we note that the direction of 
the axis of the plume is not consistent with the directions of the normals to the 
surface, but is offset by ~27°. It is, however, in much better accord with the 
direction of the normal to the long axis that passes through the sub-solar point 
(dotted line). This direction would be the direction of gravity if the gravitational 
potential is flat along a meridian as suggested in Thomas et al. (2012). The rectangle 
denotes the footprint of the IR spectrometer pixel. Note that the orientation of the 
nucleus in these images is not consistent with other images shown in this paper, 
which are orientated as in the original camera frame. In this figure the nucleus has 
been rotated anti-clockwise by 90° to match the convention used in A’Hearn et al. 
(2011) that puts the direction of the Sun approximately to the right. 

sion of water ice and C0 2 from this region). It is beyond the scope 
of this paper to discuss the details of processes in the inner coma, 
nevertheless we would note that there could be a temperature ef- 
fect at work here with the sub-solar region at the waist being con- 
siderably hotter than the region at the small end of the nucleus, 
which may be cooled by localized vaporization of C0 2 . Some indi- 
cation of this cooling is seen in the effective temperature maps of 
the nucleus in Groussin et al. (2012). 
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